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Abstract — Polymer flooding can not only accelerate the rate 
of oil production but also improve the efficiency of oil 
production. The finite volume method is used to study flow 
characteristics of viscoelastic polymer in microchannel. The 
results showed that: the displacement efficiency of the 
viscoelastic polymer is better than that of water. The greater the 
viscosity of the polymer, the greater the drag force rid of 
residual oil and the displacement efficiency increases as well. As 
the polymer’s viscosity increasing, the flow region of concave 
corner may increase but the retention area of residual oil may 
be decreased, the spread area of the polymer in the voids 
increase and the displacement efficiency increases. The results 
can provide a theoretical basis for oilfield polymer selection and 
optimization. 

Index Terms — Finite volume method, Numerical 

simulation, Viscoelastic properties, MicroChannel 

I. Introduction 

Oil is a nonrenewable resource and an important economic 
lifeline of the country's economic development. At present, 
most oilfields of the world has entered the later stages of 
water flooding development, which lead to water content 
increased and the production declined, capital investment and 
costs increased, economic indicators and technical indicators 
increased as well. Most of the remaining oil and residual oil in 
the reservoir layer could not be exploited due to differences 
between oil and water in nature, uneven distribution of oil 
layer and other factors. Water flooding is very difficult to 
open up the remaining oil and residual oil in the high 
water-cut stage. How to exploit in the cost-effective method 
is a topic which is unremitting research of the national 
reservoir engineering specialists. After years of research 
work, scientists have proposed polymer flooding, which is 
not only different from water flooding, but also can improve 
the displacement efficiency of flooding technology. It is 
applied to actual industrial production, that has achieved 
some good results [1-3]. Tang Guoqing [4] observed and 
studied the polymer flooding mechanism by using the 
microchannel model, the results showed that: the polymer 
solution can expand swept volume and control the flow ratio 
through streaming around and control profile .In addition, the 
effect of adsorption and viscous resistance and the increased 
driving pressure can improve the displacement efficiency. 
Xia Huifen et al [5-6] obtained the mechanism of improving 
microchannel displacement efficiency of polymer solution by 
studying the experimental that the viscoelastic polymer 
solution can not only reduce the water flooding residual oil, 
but also can pull residual oil out from microchannel. In 
addition, the part of the residual oil was pulled into oil slick 
which forms a new oil flow by the polymer solution. Then the 


discontinuous residual oil was gathered a flowable oil by the 
normal force of the polymer solution and ultimately the 
microscopic oil displacement efficiency was improved. This 
paper study the flow characteristics of viscoelastic polymer in 
the micro tract with the finite element volume method and 
obtain its characteristic , which is significant for further 
deepen the understanding and rich oil displacement 
mechanism and guide oilfield polymer flooding large-scale 
use and enhance recovery factor. 

II. Problem description of viscoelastic polymer flow 

The real pore shape in the residual oil reservoir is 
different and polymer is viscoelastic, which makes numerical 
calculation complex and solve difficult, the researcher 
simplify process of microscopic voids In order to facilitate 
the calculation of study, microscopic voids are simplified to 
some easier models ,for example, sudden expansion 
symmetric model, sudden contraction symmetric model , 
micro symmetric model with a blind end and so on[7]. The 
micro channel flow problems viscoelastic polymer in both 
sudden contraction sudden enlargement stage another stage is 
studied in this paper, as shown in Figure l.the ratio between 
the width of the constriction and the width of the main 
channel is 1:4, length ratio of 1: 1[8]. Take A = 10pm, B = 
2.5pm as example. The convex angular position corresponds 
to dead oil area of microscopic pores, including both shear 
flow and extensional flow, which is the focus location of this 
study. The expansion - contraction - dilation tube model is 
characterized by simple geometry of the fluid flow area, 
which includes a variety of process flow rheological behavior 
of the flow process. The flow of the polymer solution is a 
simple shear flow on away from the upstream and 
downstream of die wall of the contraction zone. The fluid 
presents a pure extensional flow in the center line of the 
contraction zone. The flow includes both shear flow and 
stretched flow in lobes location, which can fully reflect the 
viscoelastic properties of the polymer solution in the micro 
mutation interface [9-11]. 
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Fig. 1 -flow channel model 


This paper select the Bird-Carreau constitutive model to 
describe the rheological properties of viscoelastic polymer 
with the first normal stress difference as the main feature in 
the reservoir conditions. This constitutive model , which 
considered viscosity and elasticity of the polymer solution, 
can accurately reflect the fundamental rheological properties 
of the polymer solution. Under actual flow conditions, the 
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polymer solution flow in the pore is isothermal, 
incompressible and unsteady flow. Excluding the mass force, 
the control equation is as follows: 


are 0.0008 pa-s, 0.236s, 0.4627s 1 . Fig 2 shows the streamline 
diagram of viscoelastic polymer in microchannel. In the same 
condition, fig 3 shows the streamline diagram of water. 


1) Continuity equations: 
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Where 

T ; . = Stress tensor, r ' r were the stress of the first 

1 1 aa yy 

law and the stress of the second law and T xy was the Shearing 
stress, pa; 

X = relaxation time , s 

7J = ju the apparent viscosity of the fluid , m pa- s 
U - velocity, m/s 

U x = the X direction of the velocity, m/s 
U y = the Y direction of the velocity, m/s 
p = fluid density, kg/m 
p = pressure, pa 
T = temperature , T 
k = Heat transfer coefficient of the fluid 
c p = Specific heat capacity. 

2) Bird-Carreau Constitutive model equation: 

v n-\ 

+ A 2 r 2 )~ 

Where 

7j 0 = Zero shear viscosity , pa- s ; 

Tj^ = Infinite shear viscosity , pa-s. 
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Fig. 2-streamline contour diagram of viscoelastic polymer 
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Fig. 3-streamline contour diagram of the water 
Fig 2 and Fig 3 show the streamlines of viscoelastic 
polymer solution in inlet contractive stage is similar to that of 
water. A cone border is formed due to the sudden contraction 
of flow channel interface, which makes the flow lines no 
longer parallel. Two whirlpools of certain rules generates on 
the corner, whereas the spread area of viscoelastic polymer is 
larger than that of water. The streamlines deformed again 
because the channel interface is enlarged, when the polymer 
solution extrudes the throat, then enters into the large pores 
through the equal sectional holes of flow channel, and the 
extruded diameter of viscoelastic polymer solution is greater 
than that of the water in varying degrees. Fig 2 and fig 3 also 
show that water generates significant whirlpool around the 
corner, however, viscoelastic polymer solution doesn’t 
generate whirlpool. The reason of the phenomenon is 
swelling effect extrusion of viscoelastic polymer solution or 
Buras effect [12]. 


Viscoelastic polymer has some features that are storage 
modulus, relaxation time, the first normal stress difference, 
Wesson Borg Number (We = 2v/Z), Deborah Number (De = 
2 / 1) and so on. In this paper, polymer viscoelastic is only 
characterized by relaxation time and Wesson Borg number. 

The viscosity of viscoelastic polymer ranges from 0.0008 
pa-s to 0.0543 pa-s, The density of the flowing medium is 
968.8kg/m ,the density of water is 990.1kg/m ,the viscosity 
of water is 0.623 mpa-s. Inlet is defined as the flow which is 
8.0xl0-6m /s and outlet is defined as the outflow. The model 
was meshed and the boundary conditions were set. The flow 
properties of the polymer are described by the Bird-Carreau 
constitutive model and the fluid control equations, and the 
model is solved by finite volume method. 

III. Flow characteristics of viscoelastic polymer in 

MICROCHANNEL ANALYSIS 
A. Effect of Viscosity 

The viscosity of viscoelastic polymer, the relaxation time 
and the shear rate are considered in this study, and the values 


3.0 ■ 

▲ 

at 

2.5 - " 

~ 2.0 - / A 

2. A 

; 1 5 - 

3 / 

0 ) A 

£ 1.0 - A 

Cl A 

0 . 5 - 

0.0 - A 

I I I I I I I I I I I I I 

0.00 0.01 0.02 0.03 0.04 0.05 0 06 
viscosity(pas) 

Fig.4-pressure change graph with the change of viscosity 
Fig 4 shows that the flowing pressure of the polymer 
solution in the microchannel changes with the change of 
viscosity. The polymer pressure increases as the viscosity 
increasing. Fong-term heoretical study and experimental 
practice have proved that the pressure different between inlet 
and outlet not only increase when the fluid viscoelastic 
increase, but also the flowing resistance also increases. The 
force that generated by the contact portion between polymer 
and wall surface increase with the viscosity increasing and the 
flow resistance increasing. When there are residual oil in the 
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wall surface, the pulling force of polymer to residual oil will 
increase in flowing process, then that lead to the residual oil 
of pores reduce. As viscosity of the polymer increasing, the 
effect that the polymer solution carry residual oil of pores gets 
better, and the displacement efficiency gets higher. The 
results accord with the research [13] that Xia Huifen studies 
the influence of viscoelastic polymer solution to oil 
displacement efficiency. 


B. Effect of relaxation time 

Fig 5 shows the streamline contour diagram of viscoelastic 
polymer at different relaxation time. 
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Fig. 5 -polymer streamlines contour diagram at different 

relaxation time 

Fig 5 shows that the whirlpool gradually decreases for the 
same viscosity polymer in the inlet contraction state as the 
relaxation time increasing, when the fluid flows from micro 
channels to the corner, and the flow line gradually close to the 
big gap of the channel. As the fluid elasticity increasing, the 
fluid agitation in the corner of flow path occurred more 
severely, and formed the vortex that can transfer energy to 
residual oil on the corner and make residual oil into movable 
oil, then it is flooded out. The spread area of polymer 
flooding will increase, and the area of residual oil in the 
concave corner becomes smaller and flowable area becomes 
larger. This experiment further indicates that the oil 
displacement efficiency of viscoelastic fluid will be increase 
with the increase of fluid elasticity and spread area. The 
experiment that Xue Xinsheng[14] studied the blind-side 
displacement experiment by using microscopic displacement 
device also appeared the vortex, which confirmed each other. 
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Fig. 6-pressure difference varies with polymer relaxation time 

graph 

Fig 6 shows that the pressure of the polymer solution 
between inlet and outlet in the pores reduces, as the relaxation 
time increasing. The needed pressure of large elastic fluid in 
the same speed and viscosity stage is smaller, which indicates 
that the injectivity of large elastic polymer is better. As the 
polymer elasticity increasing, viscoelastic polymer has still 
an remaining recoverable elastic modulus , which makes the 
volume expansion of fluid become large, improving the 
spread area of the fluid in the pores, then reach the effect 
improving the efficiency. 

IV. CONCLUSIONS 

A detailed CFD model for the flow of polymer solution 
was successfully implemented with finite volume method. 
This method that improve the accuracy and convergence, is 
very suitable to analysis the flow characteristics of large 
viscoelastic polymer solution. This model can accurately 
reflect the actual condition of reservoir condition. The 
method can also be extended to more the constitute equations, 
for example, Maxwell model and Oldroyd-B model or other 
Simplified model of porous media. 

This model provides detailed information about the flow 
characteristics of viscoelastic polymer for different viscosity 
and different relaxation time. The analyzed results show that 
the higher viscosity polymer lead to higher drag force and the 
higher ability to carry the residual oil and higher oil 
displacement efficiency in the above model. And the spread 
area of polymer flooding is larger than that of water and it’s 
displacement effect is better than that of water flooding, the 
results are consistent with experiments of Xia Huifen et al 
( 2002 ) 

The agitation in the" Blind Side" generates the vortex that 
can make residual oil into the movable oil to get rid of corner, 
as elastic polymer increasing. The area of residual oil is 
reduced. Flooding area of residual oil is larger. And the 
injectivity of polymer get better, then the spread area of 
polymer in the gap model increases .The efficiency of 
displacement gets higher improve the efficiency of polymer 
flooding. The results are consistent with experiments of Xue 
Xinsheng et al (2005). 



References 

[1] P.P. Shen, C.Z. Yan and S.Y. Yuan. EOR polymer flooding [M]. 
Beijing: Petroleum Industry Publishing House,2006:l-2 

[2] Q.L. Gang. A dissertation on China recovery[J] .OGRT, 1998,5(4): 1-7 


29 


www.ijeart.com 


Flow characteristics of viscoelastic polymer in microchannel 


[3] J.H. Yang. Status and developing tendency of EOR research [J], 
Petroleum Geology and Recovery Efficiency, 1999, 6(4): 1-5. 

[4] Z.H. Zhang, Y.C. Xue, C.H. Xie and J.C. Wu. Operator-splitting 
matched artificial dispersivity method in solving three-dimensional 
advection dominant advection-dispersion problems [J]. Chinese 

journal of computational physics, 1994, ll(l):78-79. 

[5] D. Camilleri , S. Engelson, L. W. Lake, E.C. Lin, T. Ohnos, G. Pope et 
al. Description Of an improved Compositional Micellar Polymer 

Simulator. SPE reservoir engineering (Nov, 1987): 427-432 

[6] J.J Liu, D Feng, W.D Ma, Z.Q Yan and J Fu. Study on Polymer 
flooding of 1-3 Sand of sheng-sha sec in shengtuo oil fieldfJ] . Oil-gas 

Field Surface Engineering,2002;21(2):45-47. 

[7] H. J. Yin, H. M Jiang and C. Y Su. Flow behavior of viscoelastic 
polymer solution in the expansion chanel [J] .Acta Polymerica Sinica, 

2009(6): 519-523 

[8] H. Zhang, H.J. Yin and H.Y. Zhong. Finite volume method of the 
microfiltration of viscoelastic polymer solution[J] .DaQing Petroleum 
Geology & Development 

[9] E. F. Block, F. Clvan. Prorous-media momentum equation for highly 
accelerated flow. SPE reservoir engineering (Aug. 1988): 1048-1052 

[10] J. A. Deiber, W. R. Schowalter. Modeling the flow of viscoelastic 
fluids through porous media. ALCHEJ, 1981, 27(6): 912-920 

[11] R. J. Marshall, A. B. Metzner. Flow of viscoelastic fluids through 
porous media. I&EC fundamentals » 1967: 393-400 

[12] H. F. Xia. The percolation theory of Viscoelastic polymer solution 
and its application [M] .Beijing: Petroleum Industry Publishing 
House, 2002 

[13] H. F Xia, D. M Wang, G. Wang and F. S. Kong. Elastic behavior of 
polymer solution to residual oil at blind end [J]. Acta Pet Rolesi 
Sinica, 2006 

[14] X.S. Xue. Effect of molecular structure on associating polymer 
viscoelastic and its seepage flow characteristics in the porous media 
.Cheng Du; Southwest Petroleum University (2005) 

Zheng Jipeng is a Master of Southwest Petroleum University. His major 
is Power Engineering and Engineering Thermal Physics and engages in 
modern design and numerical research of Petroleum machinery. 

Ye Yulin holds a bachelor degree in mechanical engineering from 
Sichuan Agriculture University (China). Currently, he is completing his 
Master of Mechanical Engineering degree at Institute of Mechanical 
Engineering of Southwest Petroleum University. He is mainly engaged in 
modern design and numerical research of machinery. 

Yan Liuying holds a bachelor degree of Southwest University of 
Science and Technology in mechanical engineering. She is mainly engaged 
in the work of petroleum machinery design and control. 


30 


www.ijeart.com 


